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Abstract. We present an analysis of stellar mass estimates for a sample of 25000 galaxies from the COMBO-17 
survey over the interval 0.2 < z < 1.0. We have developed, implemented, and tested a new method of estimating 
stellar mass-to-light ratios, which relies on redshift and spectral energy distribution (SED) classification from 5 
broadband and 12 medium band filters. We find that the majority (> 60%) of massive galaxies with M, > 10^^ Mq 
at all 2 < 1 are non-star-forming; blue star-forming galaxies dominate at lower masses. We have used these mass 
estimates to explore the evolution of the stellar mass function since z = 1. We find that the total stellar mass 
density of the universe has roughly doubled since z ~ 1. Our measurements are consistent with other measurements 
of the growth of stellar mass with cosmic time and with estimates of the time evolution of the cosmic star formation 
rate. Intriguingly, the integrated stellar mass of blue galaxies with young stars has not significantly changed since 
2^1, even though these galaxies host the majority of the star formation: instead, the growth of the total stellar 
mass density is dominated by the growth of the total mass in the largely passive galaxies on the red sequence. 
Key words. Galaxies: mass function - Galaxies: evolution 
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c/3 1. Introduction In this context, an empirical assessment of the build- 
J^ ' up of stellar mass from early epochs to the present is an 
^ ' Current models for the formation and evolution of cos- important constraint. To explore the efficiency of star- 
■^ mic structure and galaxies provide us with a picture of formation and feed-back on different mass scales, one 
rN hierarchical structure formation, where large objects are needs to derive not only the integrated stellar mass density 
C^ built-up successively from the merging of smaller, previ- {p*){z) as a function of epoch, but also the evolution of 
' ously formed objects (e.g. Blumenthal et al 1984; Davis the stellar mass function of galaxies. At the present day, it 
et al. 1985) and from the more gradual inflow of gas (e.g. has become clear that stellar mass correlates strongly with 
White and Rees, 1978; Murali et al. 2002). The dynam- a broad range of galaxy properties, such as mass density, 
ical aspects of CDM-based galaxy formation can by now mean stellar population age, metallicity, or alpha enhance- 
be well described either by numerical ab-initio simulations ment (e.g.. Bender, Burstein & Faber 1993; Kauffmann et 
(e.g. Springel & Hernquist 2003; Weinberg et al. 2004) or al. 2003; Tremonti et al. 2004; Thomas et al. 2005). Study 
by semi-analytic modehng (Kauffmann et al. 1993; Cole of how some of these scaling relations evolve with redshift 
et al. 1994; Cole et al. 2000; Somerville et al. 2001). Yet, wiU give still keener insight into star formation, feedback, 
when, how and how efficiently stars formed in the course and galaxy assembly through merging, 
of this evolution appears to depend sensitively on regu- 
lating feed-back and on sub-grid (i.e. unresolved) physical Over the last decade (e.g., Lilly et al. 1996; Madau 
processes; consequently, the current generation of models et al. 1996) the global build-up of stellar mass has been 
has limited power to robustly address this most central traced by estimating the mean star-formation rate (SFR) 
question of galaxy formation. as a function of cosmic epoch, (SFR(z)). These estimates 

have been based on either the UV luminosity (e.g. Steidel 
et al. 1999; Giavalisco et al. 2004, Schiminovich ct al. 



Send offprint requests to: K. Meisenheimer, e-mail: 2005), on the resulting emission hue excitation (e.g. Yan 
meisenheimer@mpia.de et al. 1999; Hippelein et al. 2003; Brinchmann et al. 2004) 



Borch et al.: The stellar masses of 25000 galaxies at 0.2 < z < 1.0 



or on thermal infrared or radio emission (e.g., Flores et al. 
1999; Haarsma et al. 2000; Le Floc'h et al. 2005). The UV- 
and line emission-derived estimates require large and un- 
certain dust extinction corrections (e.g., Meurer et al.l999; 
Adelberger and Steidel 2000); thermal infrared and radio 
emission are less dust-sensitive but are currently limited 
in sensitivity to probing intensely star-forming galaxies at 
intermediate and high redshift. The latter approach there- 
fore requires large corrections for the faint end of the IR 
or radio luminosity function to estimate (SFR(z)). Both 
approaches rely quite sensitively on a prior assumption 
for the stellar initial mass function, as the bolometric lu- 
minosity is dominated by the UV radiation from ^ SM© 
stars, while the mass budget is made up by stars of so- 
lar mass and below. Finally, even perfect knowledge of 
the star formation density function at different redshifts, 
/9(SFR, z), does not contain enough information to predict 
the observable stellar mass function of galaxies, p(M*, z). 

A conceptually obvious, but practically quite difficult, 
path to bypass these problems in estimating the build-up 
of stellar mass, is to estimate p(M*, z) directly from stel- 
lar mass measurements at different redshifts. Two obser- 
vational avenues present themselves, dynamical mass esti- 
mates from kinematic measurements at different redshifts, 
and estimates of the stellar mass-to-light ratio {M / L) 
based on the spectral energy distribution (SED) of a stel- 
lar population. Dynamical mass estimates from molecular 
gas kinematics have now been demonstrated to redshifts 
z > 2 (e.g. Neri et al. 2003), but will remain restricted 
in the immediate future to exploring the relatively small 
fraction of massive, gas-rich galaxies. Kinematics of galax- 
ies to 2; > 1 are now routinely being measured, both from 
absorption lines (e.g., Holden et al. 2005; van der Wei et al. 
2005) and emission line rotation curves (Vogt et al. 1996, 
1997; Bohm et al. 2004; Consehce et al. 2005); current 
published sample sizes are in the 10s to 100s. 

Stellar luminosity is the simplest observational proxy 
for stellar mass, but cannot offer stellar masses with rele- 
vant precisions (better than about a factor of two) , as the 
rest-frame optical M / L values of stellar populations with 
different SFHs can vary by more than an order of magni- 
tude. Rest-frame near-infrared M / L ratios vary less as a 
function of star formation history; nevertheless, different 
plausible star formation histories differ still substantially 
in their rest-frame infrared M / L (by a factor of four or 
more). 

However, the luminous properties of galaxies can be 
used to estimate stellar masses. Low-resolution (R ^ 
3 — 30) galaxy spectral energy distributions (SEDs) do 
not contain enough information to robustly infer the SFH, 
metallicity, or dust content of a galaxy. Older stellar ages, 
higher metallicity, or more dust extinction all make SEDs 
redder. Yet, all these effects not only redden the SED, 
but also lower the emergent optical flux (at a given mass) 
by similar amounts. Hence, stellar M/ L can be estimated 
from SED-template fits, the strength of the 4000A- break 
(D4000) or even the (B-V)rost color, for a wide range of 
the above parameters (e.g., Brinchmann & Ellis 2000; Bell 



& de Jong 2001; Papovich et al. 2001; Kauffmann et al. 
2003) . If spectra are available, this approach can be refined 
to include two or more spectral diagnostics, e.g. D4000 
and EW(H(5) (Kauffmann et al 2003); indeed, recent works 
make use of higher-resolution spectra to robustly estimate 
stellar mass (e.g., Panter et al. 2004). Stellar masses esti- 
mated in these ways are consistent at the ~ 0.1 dex level 
(e.g., Bell et al. 2003; Drory et al. 2004b; Salim et al. 2005) 
and correlate strongly with dynamically-measured masses 
(e.g., Bell & de Jong 2001; Drory et al. 2004a; van der 
Wei et al. 2005; Cappellari et al. 2005). One has to be 
aware, however, that very dusty galaxies or galaxies with 
SFHs featuring by large and recent (less than 2 Gyr ago) 
bursts of star formation will have poorly-measured masses 
(^ 0.3 dex error) if not derived from template spectra ac- 
counting explicitly for these effects. 

SED or spectra-based M/L estimates have been used 
to describe the local galaxy population (Cole et al. 2001; 
Bell et al. 2003; Kauffmann et al. 2003), based on samples 
of many thousands of galaxies, and on modest samples 
of distant galaxies z > 1, where deep near-IR data have 
permitted measurement of the optical rest-frame SEDs 
(e.g. Papovich et al. 2001, Rudnick et al. 2003, Drory 
et al. 2004a, 2005; Fontana et al. 2004). Studies of dis- 
tant z > 0.5 samples have been hampered by a number 
of largely inevitable difficulties: small number statistics 
for the most massive galaxies, field-to-field variations (as 
quantified elegantly by Drory et al. 2004a and Somerville 
et al. 2004), and increased stellar mass uncertainties due 
to limited spectral coverage and/or the use of photometric 
redshifts. 

The purpose of the present paper is to derive SED- 
bascd stellar mass estimates for all galaxies in the 
COMBO-17 sample (Wolf et al. 2001,2003). We quantify 
field-to- field, number statistics, and photometric redshift 
uncertainties as carefully as possible. We then use those 
mass estimates for a comprehensive estimate of the global 
stellar mass density in the redshift range 1 > z > 0.2, 
as well as of the galaxy mass function over the same pe- 
riod. This paper is organized as follows. Section [5] briefly 
recapitulates the data base of the COMBO-17 survey and 
outlines its redshift and SED determination. Section|31de- 
scribes the methodology for estimating galaxy masses and 
discusses the errors involved. Section^presents our results 
on the evolution of the stellar mass in galaxies since z = 1. 
In Section we discuss these results in respect to previ- 
ous work, the cosmic star formation history and models of 
galaxy formation. We assume Hq — 70/i7okms^^ Mpc~^; 
flA = 0.7 and flm = 0.3 and adopt hro = 1. 

2. The galaxy sample 

Our analysis is based on the optical multi-color images in 
three field s (CDFS, Abe ll 901, llh field) of the COMBO- 
17 survey IjWolf et al. 2003:) . We use the photometric cat- 
alog of Wolf et al. (2003), in the revised version presented 
by Wolf et al. (2004) to which we refer for further details. 
Each square field covers ^ 950 D' on the sky, correspond- 
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ing to a co-moving transverse size of ~ 22 Mpc at z~ 0.7. 
The three fields are so widely separated as to be com- 
pletely uncorrelated. The combination of deep observa- 
tions in 5 broad bands (UBVRI) and 12 medium band fil- 
ters between 400 and 930 nm, allows us to determine very 
accurate photometric redshifts (see Section im below and 
Wolf et al. 2004) and accurate UV-optical spectral energy 
distributions of galaxies fSection l2.2|l . Together, the three 
fields contain more than 26 000 galaxies with redshift esti- 
mates out io z = 1.07 which is the limit of reliable, precise 
redshift estimates based on optical filters only. 

The multi-color photometry can be used to a) clas- 
sify each object as star, galaxy, QSO or indeterminate, 
h) estimate its redshift, e.g. if it is a galaxy, and c) fit a 
template SED to the fiuxpoints, and derive physical pa- 
rameters, such as the stellar M/L. In principle, such a 
process can be carried out in one step (e.g., Rudnick et al. 
2003) where classification, redshift and stellar mass are es- 
timated jointly. However, experimentation with COMBO- 
17 data showed that physically- motivated template se- 
quences with composite (old-|-young) stellar populations 
gave poorer redshift accuracy than the Wolf et al. (2004) 
dust-reddened single-burst template set. Consequently, in 
this paper we have adopted a two step classification pro- 
cess. First, we classify the objects and estimate redshifts 
using dust-reddened single burst templates (Wolf et al. 
2004), which have been proven to give the best photomet- 
ric redshifts as tested against a large spectroscopic train- 
ing dataset, but clearly cannot provide the most realistic 
stellar masses. Therefore, the stellar masses of all galaxies 
and their uncertainties are obtained by fitting the opti- 
cal 17-band SEDs with a new set of templates with more 
plausible SFHs. 

2.1. Object classification and redshift estimates 

In COMBO-17, objects are classified by their location in 
the 16-dimensional color space IjWolf et al. 200 l|l . We use 
four broad-band color indices: U — B, B — V, V — R, 
R — I, and 12 medium band indices which always relate to 
the closest broad band filter. Color libraries of four object 
classes, namely main sequence stars, white dwarfs, galax- 
ies and quasars, are used as templates. The galaxy tem- 
plates span the redshift range < z < 1.4, the quasars 
the range < z < 5. In the first step, each object is 
assigned to one of these four types by comparing its ob- 
served colors (and their errors) to the color distribution 
of the templates. For an unambiguous classification, we 
require a relative probability oi pi/J^Pi ^ 0.75. In fact, 
98% of all objects at R < 24 are unequivocally assigned to 
one class by this requirement. The object classification re- 
lies primarily on colors; morphological information is used 
only for objects with Pi/J^Pi < ^-"^^ ^^^ ^^ classes. 

For galaxies and quasars, we derive the joint proba- 
bility of a given redshift and a given rest-frame spectral 
energy distribution (SED), accounting for photometric er- 
rors. The galaxy SEDs are represented by a set of simple 
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Fig. 1. Parameterization of the star formation history in 
our templates. An exponentially decaying star formation 
rate simulates an initial burst (component a). A second 
component b provides some constant star formation since 
the initial burst. For starburst galaxies a second burst in 
the recent past is assumed (component c). 

templates (single-burst stellar populations with redden- 
ing) that was tuned to provide the most robust redshift 
estimates through comparison with spectroscopic redshifts 
(see IWolf et al. 2004|l . This procedure provides a mini- 
mum error variance estimation of both the redshift and 
the SED template. For this paper, galaxies with well- 
measured minimum-error- variance redshifts are used, with 
an additional R-band aperture magnitude cut Rap < 24 to 
ensure that only galaxies with reasonably reliable redshifts 
are included in the analysis. 

Wolf et al. (2004) have shown that in the redshift range 
O.K z < 1.05 the redshift accuracy of COMBO-17 is weU 
approximated by 



a, = 0.007(1 + z)\/l + lOO-8(R-2i.6)^ 

where R is the central R-band magnitude of the galaxy 
(measured in a 1'.'5 aperture), and does only weakly de- 
pend on galaxy color (SED). 

2.2. Determining tiie SED of the galaxies 

Obviously, the single burst templates used for the red- 
shift estimates hardly can account for the wide variety of 
SFHs in galaxies. Detailed observations of resolved stel- 
lar populations (e.g., Rocha-Pinto et al. 2000; Harris & 
Zaritsky 2004) and integrated spectra (e.g., Trager et al. 
2000; Proctor & Sansom 2002; Panter et al. 2004) show 
that most galaxies have formed stars over extended pe- 
riods — oftentimes over almost a Hubble time. As old 
stars usually dominate the mass while young stars often 
dominate the optical flux, it is essential for estimating the 
mass-to- light ratio M/L to constrain their relative amount 
as accurately as possible from the observed SED. 

Adopting the classification and redshifts from the 
simple templates, we estimate (M/L)^ with a second 
SED template set that reflects more plausible SFHs. 
Specifically, we have devised a library of 100 spec- 
tral energy distributions built with the PEGASE-code 
p'ioc &: Rocca-Volmerange 1997| ), sorted from red to blue 
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by an SED index 5, < 5 < 99, which reproduces the 
sequence of mean UV-optical spectra of nearby galaxies 
collected by Kinney et al. (1996) ^ The underlying SFHs 



have been parameterized by the three-component model 
depicted in Figure ^ For the " old" population, we as- 
sume an exponentially decaying SFH (with e-folding time 
1 Gyr) several billion years ago (component a), plus a 
constant star formation rate since the initial burst (com- 
ponent h). To reproduce the empirical SEDs between E 
galaxies and Sb galaxies in Kinney et al. (1996) we con- 



tinuously decrease the age of initial burst by a factor 
of 2.5 and increase the h component to a maximum of 
h/a = 0.16^. For even bluer galaxies (the "starburst galax- 
ies" SB6 to SBl in Kinney et al. (1996) I we include a sec- 
ond burst in the recent past (60 million years ago, com- 
ponent c). Its exponential decay time is kept constant, 
Tc = 10* yrs. Bluer galaxies are generated by increasing 
the relative amount of recent star formation c/a according 
to the prescription 



c/a ^ 0.002 e 



S/17.2 



0.021, 



in the range 40 < S* < 99. 

For component (a) we chose to follow the chemical evo- 
lution of the model galaxy using a closed box model, start- 
ing at 1/20 solar metallicity, as implemented in PEGASE. 
The metallicity of components (b) and (c) are chosen to 
be constant at solar metallicity. As we assume no infall 
and no galactic winds and start at 1/20 solar metallic- 
ity, our models will exhibit a significant G-dwarf prob- 
lem. Furthermore, super-solar metallicities cannot be eas- 
ily reached by this model; thus, the reddest SEDs (E, SO 
type) can only be reproduced using unreasonably old ages 
of 20 Gyrs. Therefore, we regard the "ages" as meaningful 
only in a relative sense, with the age of the component a of 
the Sb galaxy found to be about QAxtmax = 8 Gyr (where 
tmax refers to the age of the oldest stars in E/SO galaxies)^. 
We have adopted a Kroupa IMF ( |Kroupa et al. 1993| ) and 
the mass regime 0.1-120 Mq. The stellar M/L values 
and colors from such an IMF are very similar to those 
of a Kroupa (2001) IMF or a Chabrier (2003) IMF. Dust 
extinction is crudely accounted for using PEGASE (see 
Fioc & Rocca-Volmerange 1997 for details). In the code 
the dust optical depth is estimated from the gas mass 



^ In this way, the new library closely matches the template 
library which has been used in early implementations of the 
multi-color classification (see Wolf et al. 2001; 2003). 

^ Parameters a, fe, and c are defined in terms of total stellar 
mass contained in each component; see, e.g.. Equations 4 and 
5. 

^ It is important to note that, owing to the age/metallicity 
degneracy for old stellar populations, the modest changes in 
stellar M/L caused by the old ages of the S < 40 templates 
are cancelled out by the substantially sub-solar metallicities 
of these templates. Thus, the template color-stellar M/L ratio 
correlation is in fact very similar to a solar metallicity, 7 Gyr old 
population for the 5* < 40 templates; such ages and metallici- 
ties are reasonably appropriate for modelling the 0.2 < z < 1.0 
galaxy population. 
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Fig. 2. An example for the template spectra for a star- 
burst galaxy. The original Kinney et al. spectrum (SB-2) 
is shown and our best fit from the PEGASE spectra. It is 
composed of an old stellar population {a -\- h component 
in Fig.^ and a young stellar population (c component). 

and metallicity, scattering is accounted for crudely fol- 
lowing Calzetti et al. (1994). For component a we choose 
spheroidal geometry, while a plane-parallel slab geome- 
try, averaged over all inclinations, is chosen for compo- 
nents b and c. This dust prescription reproduces many of 
the broad trends observed in local galaxies (e.g., Wang & 
Heckman 1996) reasonably well**. 

In order to be suitable for multi-color classification, the 
one-dimensional library S{{) . . . 99) = 5(a, 5, c) has to be 
monotonic both in color space and in the space of its astro- 
physical parameters: b/a,c/a and the age of the old stellar 
population. Figure |2 shows an example of one starburst 
galaxy from Kinney et al. (1996) and its composition from 
an old plus a young stellar population. For comparison 
with previous work we characterize the SED types S also 
by their rest-frame colors U-V and B-V (see Fig.O. 

The best-fitting SED is obtained by comparing the ob- 
served colors of each galaxy with the colors of the SED 
library spectra, shifted to the most likely redshift z as 
measured using the Wolf et al. (2004) classifier (Section 
EH). 

3. Estimating stellar masses 

Based on the (M/L)* of the best-fitting new SEDs, stel- 
lar masses can now be estimated. We choose to normal- 
ize the stellar mass by the reddest well-measured filter to 
minimize the influence of bursts of star formation. The 
medium-band filter which provides the best compromise 



^ The color effects of age, metallicity and dust are largely de- 
generate when estimating stellar M/L from broad-band SEDs, 
therefore the detailed choice of this or a different dust model 
does not significantly affect our derived stellar masses. 
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Fig. 3. Rest-frame colors U-V (solid line) and B-V (dotted 
line) as function of the SED type S in our library of galaxy 
templates. 

between longest observed wavelength and signal-to-noise 
ratio is 816/21, with a central wavelength of 816 nm and 
a width of 21 nm. The 816/21 filter has a limiting mag- 
nitude of 22.9 (lOcr) in all of our fields. The lower edge of 
the 816/21 filter, Amin = 804 nm, sets the maximum red- 
shift where the stellar mass can be measured reliably: for 
a significant light contribution from the old stellar popu- 
lation, one needs to make sure that only flux at rest-frame 
wavelengths Arcst > 400 nm is sampled. Hence, we restrict 
our current investigation to redshifts z < 1.0. 

In practice, we convert the observed quantities 
z, S, /gie of each individual galaxy into stellar mass in 
the following way: 

(1) The luminosity delivered from the PEGASE code, 
which is normalized to 1 Mq, is used to calculate the ap- 
parent flux fx. For simplicity, we first consider only one 
component (luminosity Fx). In this case the observed flux 
of a galaxy of stellar mass M is given by: 



fx 



jFx-QxdX 
jQxdX 



1 M 

m Mq 



i7TDl{z){l 



2)/^7o' 



(1) 



Here m ■ Mq denotes the stellar mass of the PEGASE 
template. The integral runs from to oo. Q\ = q{X)t{X) 
is the spectral response of the system, in our case given 
by the transmission of the 816/21 filter and the quantum 
efficiency of the CCD. The factor (1 -I- z) accounts for the 
different wavelength scales of the observed flux fx and the 
rest- frame luminosity Fx- Both Fx and fx are given in 
Wm~^nm~^. 

In practice, COMBO-17 measures the central flux 
f'y,ap of each object in photon units 7m~^s~^nm~^ 
jWolf, PhD thesis| ) after convolving all images to a com- 
mon effective PSF of 1.5 arcseconds. Thus, for well re- 
solved galaxies, the central flux f-y.ap has to be corrected 
to the total flux f-y^tot by a correction factor Tjtot- It is de- 
rived for each galaxy from the SExtractor MAG-BEST mag- 
nitudes (Kron magnitudes) on our very deep R-band im- 
age: T]tot = 10'^*(^''p^^^'°"). Since we ignore wavelength 



dependent aperture effects, the apparent flux fx in equ. 
(0 is approximated by: 



/a 



he 



-/• 



he 



■^ ,tot 



~ J-jjap ' Vtot 



(2) 



With equation |(2J), setting f^^ap — fsw, and solving for- 
mula CQl for the mass one gets: 



1 M _fsi6Vtof'inDl{z)-{l + z)!ffQxdX 



tuMq 



h^,o-jFx-QxdX 



(3) 



(2) As our templates are composed of three components 
a,b,e, both the luminosity and the stellar mass have to be 
synthesized: For each SED type Si = S{a, b, c) the total 
luminosity (per Mq) is 



Fx = a- Fx,a + b ■ Fx.b + c ■ Fx., 



(4) 



Similarly, the total stellar mass results from a linear com- 
bination of the individual components with the parame- 
ters a, b and c: 



m — a ■ nia + b ■ nib + c • rric. 



(5) 



Here ma, mi, and rric are the stellar mass fractions of the 
respective components. Inserting (jSJ and Q into equation 
©, yields the final result: 



h 



M 



fsw mot ■ ^TTDlJz) .{l + z)!fjQxdX 

J2k=a,b,c TKZ I P^-k ■ Q\ dX 



(6) 



Equation @ delivers the stellar mass of any galaxy with 
observed values /sie, '7tot, for which the classification has 
succeeded in determining the redshift z and the spectral 
type S — S{a,b,c). The < 2% of galaxies for which no 
reliable redshift can be determined are faint and blue: mis- 
placing them does only affect the mass function of blue 
galaxies below 3 10^ Mq (see Fig.EjI. 

3.1. Intrinsic errors 

As obvious from equation {^, the intrinsic errors of the 
stellar masses depend both on the errors in redshift z {via 
(1 + z)D\{z)) and in SED type S {i.e. the M/L ratio: 
M/Fx). Since a red SED at lower redshift can produce 
similar colors as a blue SED at slightly higher redshift, 
the errors in z and S are correlated. To account properly 
for this dependence we calculate the covariance matrix 
for each galaxy. To this end, besides determining the SED 
type S{zq) for the optimum redshift (according section 
12.11 with its error az)j we let z vary in the interval [zq — 
2crz,zo -f 2az]. From the corresponding points Si{zi). we 
calculate the covariance matrix 



pf(^r^^' 



Pi{S,-S{z^){^-^Y^) 



C = 



Pf{J-^Y^){S,-S{zo)) P^{S,^S{zo))^ 
where each point is weighted by the probability 



P, 



exp 



1 fz,, 
2 



20 



(7) 
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This matrix has to be diagonalized to obtain independent 
new parameters C, and /i, calculated for each galaxy from 
the parameters z and S. The mass error is then estimated 
by gaussian error propagation from these new parameters 
and the error of the observed flux, ctj : 



fOM 



(DM 



/ dM 



-c^/ 



(8) 



3.2. Calibration of mass scale and systematic errors 

To investigate the global (average) properties of hundreds 
or thousands of galaxies within a redshift range, the sta- 
tistical mass error of an individual galaxy is less important 
than any systematic error in the mass scale. We identify 
three main sources of systematic errors: 

1. Overall calibration of the mass-to-light ratio scale. 

2. SED dependent errors caused by an incorrect conver- 
sion from observed colors of a galaxy into an SED S 
(i.e. essentially: assuming an incorrect SFH which like- 
wise fits the colors). 

3. Redshift dependent errors, which could be either 
caused by sampling different rest-frame wavelengths 
or by using a common set of template spectra for a 
galaxy population which evolves with cosmic time. 

In order to estimate how much these systematic effects 
might influence our results, we first compare (Fig.^ the 
color-M/L relation from Bell et al. (2003) which is an 
update of Bell & de Jong (2001), scaled to a Kroupa et 
al. (1993) IMF, with the result of our mass estimation 
for all galaxies at 0.1 < z < 0.6 (where the rest-frame 
V is sampled by the COMBO-17 photometry), and with 
i?-band aperture magnitudes Rap < 23 (to avoid exces- 
sive statistical errors) . The discrepancies both at the blue 
and red end can be easily understood by the fact, that 
our template age-metallicity combination (red: metal poor 
and unreasonably old; blue: solar metallicity and young 
^ 0.1 Gyr) should result in stellar M/L ratios at a given 
optical color which are very slightly smaller than those 
of Bell & de Jong (2001), who assumed 12 Gyr old popu- 
lations and relatively high metallicities. Notwithstanding 
this slight difference, we conclude that the stellar M/L 
ratios used in this paper are very similar to those used 
for local studies, arguing that comparisons of this work to 
z ^ studies should be reasonably robust. 

Since our stellar masses are derived from observed- 
frame optical data alone, an independent second check of 
the robustness of the stellar mass estimates as a function 
of SED type is to compare the optical-only stellar masses 
with stellar masses derived from observed-frame near- 
infrared data. We have used deep, public, fully-reduced J, 
H, and K-hand data on a 50 square arcminute sub-region 
of the CDFS taken with ISAA C at the VLT as a part of th e 
ESO Imaging Survey EIS Patziminaoglou et al.2002| ). 
Their limiting magnitudes are J< 24.8, H< 23.4 and 
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Fig. 4. V-Band mass to light ratio versus (B-V)rcst color 
of COMBO-17 galaxies. For comparison the Bell et al. 
(2003) color-M/L correlation is shown as straight line. 
The transition between early type {S < 40) and blue 
galaxies occurs at {B — l^)rest = 0.8. Also overplotted is 
the color track of a sequence of exponential star forma- 
tion histories with solar metallicity and age 7 Gyr (solid 
curve). 

Kg < 22.2 (5cr, Vega system). Although this is only a 
small fraction of one COMBO-17 field, we find a sample 
of about 400 galaxies for which both a rehable COMBO-17 
classification and near-infrared photometry are available. 

We feed the ISAAC images directly in our COMBO-17 
pipeline and re-calibrate the measured flux with respect 
to main sequence stars in the fleld.^ In order to allow 
matched photometry, the field distortions of the ISAAC 
camera have been corrected radial symmetrically, based 
on relative positions of bright objects on our deep R-band 
image and the individual J, H, K pointings. 

For the comparison, we do not re-calculate best-fit red- 
shifts and SED types using the optical and near-infrared 
data. Instead, we use the optically-derived redshifts and 
SED fits to blindly predict stellar M/L values for the 
observed-frame J, H, and K-hands, which are then mul- 
tiplied by the J, H, and A'-band photometry to derive 
new stellar masses^. That is, we do not attempt to find 



^ This re-calibration accounts only for the global zero-point 
of the EIS data. From the main sequence stars we have some 
indications that the zero-points of the 8 subfields are slightly 
inconsistent with each other. 

® Recall that the default mass estimates were derived from 
SED-derived M/Ls in the observed-frame 816/21 passband, 
nmltiplied by the 816/21 observed-frame luminosity. 



Borch et al.; The stellar masses of 25000 galaxies at 0.2 < z < 1.0 



in 



O 



0.5 







-0.5 



- 






■4-'" ++ S 








Pi.?v, 


- 










■I * ^ * *'"'"''{*'' )t*"I. +1; ■"'■*' 



o 




20 40 60 
SED 



80 



fold /frr 



Fig. 6. Inferred stellar mass-to-light ratio for starburst 
galaxies as a function of the age of the underlying old stel- 
lar population (continuous line: V-band, dashed: K-band 
M/L) . The plot show the ratio between the true M/L ra- 
tio and that inferred with our toid/tmax = 0.4 templates 
1 UU when classifying starburst galaxies for which we varied the 
age of the old stellar population. Note that both the V- 
band and K-band M/L ratio would be affected to a similar 
extent. The input template types were 50 < S" < 90. 



Fig. 5. Ratio between the stellar mass estimation based on 
the if -band and the stellar mass from the 816/21 filter as 
a function of SED type. The crosses with error bars show 
averages in SED bins (A^* — 10) and their rms variation. 



a better overall solution for fitting the template spectra 
to all observed colors (16 optical plus 3 including the 
near- infrared) . This approach has the advantage that we 
are sensitive only to the effect due to sampling different 
parts of the galaxy spectra for estimating the stellar mass. 
As the relative fraction of light from young and old stel- 
lar populations is a function of the rest-frame wavelength 
which is sampled, any inconsistencies between our synthe- 
sized spectra (and their composition) and the true spectra 
of the galaxies should become evident. 

In Figure the comparison between the mass estima- 
tions from filter 816/21 and from the i/-band is shown as 
a function of SED type. The agreement is on average much 
better than 0.1 dex (25%) for early type and spiral galaxies 
(5 < 70). Later templates — those with significant bursts 
of star formation in our template set with S* > 70 — are 
offset by up to 0.15 dex in the sense that the optically- 
derived masses overestimate the stellar mass. The object- 
by-object scatter is largely accounted for by random er- 
rors in stellar mass: more than 60% of galaxies have mass 
differences of less than the combined mass uncertainties. 
The use of J-band for the mass derivation leads to similar 
results. If X-band is used for the mass derivation, again 
galaxies at S" < 70 show consistent masses to within their 
combined uncertainties, whereas galaxies with the bluest 



stellar populations 5* > 80 show K-hand derived stellar 
masses —0.3 dex lower than the optically-derived masses. 
While this offset is significant, it is worth bearing in mind 
that the masses of the bluest, most strongly star- forming 
galaxies were expected to be systematically uncertain at 
this level (owing to the unknown frequency of bursts of 
star formation; Bell & de Jong 2001; Kauffmann et al. 
2003). Furthermore, galaxies with S" < 70 — whose stel- 
lar masses appear to be consistently derived using either 
optical or near-infrared data — form more than 80% of 
the total integrated stellar mass at all z < 1. Thus, we 
conclude that the global results presented in this paper 
would not significantly change if near-infrared data were 
used in addition to the optical data to derive the stellar 
mass function evolution. 

Although Bell & de Jong (2001) have demonstrated 
that there is an average correlation of the stellar M/L ra- 
tio and optical color for a wide variety of SFHs (see Fig. 
01 , different SFHs resulting in the same optical color might 
well differ by a factor of 2 in their M/L ratio. Our SED li- 
brary has been constructed under the assumption that all 
blue galaxies {S > 40) possess a common old stellar pop- 
ulation of age ioid = 0.4tniax- Since at z = 1 the universe 
was half its present age, this assumption might introduce 
a redshift dependent systematic error in the mass of blue 
galaxies. As a third check we have investigated this effect 
by constructing blue galaxies with various values of ioid, 
subsequently classified with our template library (constant 
told — 0.4iniax)- As One can see from Fig.El increasing the 
age of the old stellar population by a factor of 2 {i.e. from 
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about 5 to 10 Gyrs) would increase the mass in blue galax- 
ies by about 50%. On the other hand, assuming a much 
younger age, t^id = 0.2imax could lower the mass in blue 
galaxies by 40%. We conclude that this effect introduces 
a ^ 0.1 dex uncertainty in the evolution of stellar mass in 
blue galaxies in the 4 Gyr from z = 0.9 to z = 0.3 (the 
redshift range probed by our study), and a <C 0.1 dex un- 
certainty in the evolution of total stellar mass over that 
period (as blue galaxies do not dominate the total stellar 
mass at any redshift). 

4. Results 

4.1. The relationship between color and stellar mass 

FigureHshows the stellar mass estimates for the COMBO- 
17 galaxies as function of rest-frame U-V color. The solid 
black line indicates the (evolving) separation of 'red se- 
quence' from 'blue cloud' galaxies, proposed by Bell et al. 
(2004), after conversion into the M, — ([/ — F)rcst plane: 



{U - F)rcst > 0.227 logio M, - 1.16 - 0.352z, 



(9) 



Not only locally (z < 0.3) but also at high redshift z ~ 1 
the intrinsically red (and for the most part old and non 
star forming) galaxies constitute over 60% of the very mas- 
sive galaxies (M,> 10"Mq). 

4.2. Completeness limits 

It is obvious from Fig.[3 that the R-band flux limit of 
COMBO-17 leads to a mass cut (denoted roughly by the 
grey line in each panel of Fig.[7J, which does not only 
depend on redshift but quite strongly on the rest-frame 
galaxy color: the least massive red galaxies contained in 
our sample typically have 10 x higher masses than the 
least massive blue galaxies. Thus, the strict mass com- 
pleteness limit is set by the red sequence galaxies. We es- 
timate an approximate completeness limit by selecting the 
mass value above which we find 80% of our detected red 
sequence galaxies in the respective redshift bin (vertical 
dotted lines in each redshift panel of Fig.d) . These strict 
limits are very conservative; in particular, the blue galaxy 
mass function will be complete well below this limit. The 
same limits, evaluated for the upper boundary of each red- 
shift bin, are also shown in the mass function plots (see 
Figs.El&EI)- 

4.3. The evolution of the stellar mass function 

We now use the estimated stellar masses of galaxies up to 
2 = 1.0 to derive the evolution of the stellar mass func- 
tion. Since our redshift slices are thin (Az = 0.2) and the 
limiting mass is evaluated at the upper boundary of each 
bin, our samples are effectively volume limited down to 
that limiting mass and we do not use a Vmax correction. 
We apply individually-derived completeness corrections to 
each detected galaxy (see, e.g.. Wolf et al. 2003 for a full 
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Fig. 7. Stellar mass versus rest-frame U-V color in differ- 
ent redshift bins. Middle redshift of the Az = 0.1 wide 
bins is given in the upper left of each panel. The solid 
black line shows the seperation between red sequence and 
blue cloud galaxies (equation (10)). The grey line indicates 
the approximate effect of the adopted apparent magnitude 
cut {ruB, < 24) for successful redshift classification. The 
dotted lines indicate completeness limits for red sequence 
galaxies; blue cloud galaxies are complete well below this 
limit. 



description of the completeness properties of COMBO- 
17). These were derived by putting artificial galaxies into 
COMBO-17 images and applying the COMBO-17 detec- 
tion and classification pipeline. The completeness correc- 
tion was derived as a function of (z, ma.ap, U —V) by con- 
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Fig. 8. Stellar mass function of all galaxies in different redshift bins (filled circles). The error bars account for Poisson 
uncertainties and cosmic variance as estimated from the full range of field-to-field variation (dark grey) and following 
Somerville et al. (2004; light grey). The solid line represents a Schechter function fit (with constant faint end slope 
Q, assumed to be a = —1.1 following the z ~ result of Bell et al. 2003). The lower grey line shows the local mass 
function; the upper grey curve shows the local mass function convolved with a log normal distribution of width 0.3 
dex to simulate the effect of photometric redshift and stellar mass errors on the recovered mass function of high-mass 
galaxies. The vertical dotted line shows the completeness limit (for red galaxies) taken from Fig. [3 Open circles 
show comparison mass functions from Drory et al. (2004a), and the dotted line in the highest-redshift bin shows the 
Schechter fit to the stellar mass function from Fontana et al. (2004). 



structing the ratio of number of recovered galaxies to by 
the number of galaxies inserted into the same bin in red- 
shift, i?-band aperture magnitude and U — V rest-frame 
color. That is, this completeness correction corrects the 
number of detected galaxies in a given bin in redshift, 
rest-frame colour and magnitude: in practice, this correc- 
tion does not exceed 20% for objects classified as galaxies. 
It does not correct for galaxies not classified by COMBO- 
17, i.e. which are too faint to be classified or have highly 
unusual 17-passband SEDs. 

The results are shown in Fig.|Sl (total for all galax- 
ies) and Fig. El (split into blue and red galaxies). The pa- 
rameters of the Schechter fits (shown by the continuous 
and dotted lines in Figs. |H1 & EJl are given in Table ^ 
Schechter functions are fit to galaxies only above our con- 
servative completeness limit (as denoted by the dotted 
line). The data lack the depth to allow robust estimation 
of the faint end slope a. A number of strategies could have 



been adopted; we chose to assume the same values of a 
as found locally: a — —1.1 (all galaxies), a = —0.7 (red 
galaxies), and a — —1.45 (blue galaxies) following Bell et 
al. (2003)^. Mass function parameters for the local sam- 
ple are reproduced in Table ^ converted to the cosmology 
adopted in this paper and using a universal Kroupa et al. 
(1993) IMF (for reference, stellar masses calculated with 
a Chabrier 2003 IMF would be similar to within a few 
percent) . 

Figure |H1 shows clearly that the evolution of the to- 
tal mass function is relatively modest. Since the low-mass 
end is completely dominated by blue galaxies (c./. Fig.l^J, 
for which the completeness reaches to smaller masses, the 



^ Other strategies could have been adopted for constraining 
a, e.g., adopting the faint end slopes of the luminosity functions 
as derived by Willmer et al. (2005). We chose to adopt values 
of a fixed to locally-determined values to facilitate comparison 
of mass functions from z = 1 to 2 = 0. 
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Fig. 9. Stellar mass function of red (o) and blue (*) galaxies in different redshift bins. The solid black line represents 
a fit of the red galaxies with a Schechter function (with a = —0.7). The Schechter fit to the blue galaxy stellar mass 
function is shown as dotted black curve (assuming a — —1.45). The local mass functions, taken from Bell et al. (2003), 
are shown in grey. The error bars account for Poisson and ficld-to-ficld uncertainty. The vertical dotted line shows the 
completeness limit for red galaxies. 



stellar mass function continues to rise well below the for- 
mal mass limit (recall, the mass limit is defined for red 
galaxies) . We show also near- infrared derived stellar mass 
functions from MUNICS (Drory et al. 2004a; open circles) 
in the 0.4 < z < 0.6, 0.6 < z < 0.8, and 0.8 < z < 1.0 
bins, and the Schechter fit to the K20-derived mass func- 
tion of Fontana et al. (2004) in the 0.8 < z < 1.0 bin. 
Overall, the agreement between COMBO-17's and other, 
near-infrared derived z < 1 stellar mass functions is excel- 
lent if one takes into account the large error bars at the 
bright end (due to small number statistics). In a recent 
study Bundy et al. (2006) present K-band based mass es- 
timates of 8000 galaxies from the DEEP2 redshift survey 
(Davis et al. 2003). Although their split between red and 
blue galaxies, and redshift binning differ from ours, they 
find very similar results. 

4.4. The evolution of the stellar mass function of red 
and blue galaxies 

Owing to the COMBO-17's number statistics and photo- 
metric redshift quality, it is possible not only to constrain 



the total stellar mass function, but to also construct mass 
functions of the red sequence and blue cloud separately. 
The results are shown in Fig.lHl Diamonds and solid lines 
(black: fit to the data, gray: local color-split mass function) 
show the stellar mass function of red sequence galaxies. 
Asterisks and dotted lines show the stellar mass function 
of blue cloud galaxies. The vertical dotted line shows the 
position of the completeness limit for red-sequence galax- 
ies: it is clear from the data points that red sequence 
galaxies become dramatically incomplete at low masses, 
whereas the mass function of blue galaxies continues to 
rise steeply. 

Focusing on the blue cloud galaxies (asterisks and dot- 
ted lines), it is clear that there is relatively little change 
in their mass function since z ~ 1. In this framework, the 
observed rapid changes in L* in the blue galaxy luminosity 
function (e.g., Lilly et al. 1995; Wolf et al. 2003; Willmer 
et al. 2005) have to be interpreted in terms of an unchang- 
ing characteristic stellar mass M* and a rapidly evolving 
mean stellar Ad/L ratio (and therefore a rapidly evolving 
average color with redshift, as is indeed observed). 
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Table 1. Result of the Schechter fits to the mass function. 



(fy* X 10" 



logio M* 



log 10 P* 







All galaxies 




0.0 


35±4 


10.91±0.10 


-1.1 ±0.02 


8.48±0.10 


0.3 


19±9 


11.03±0.08 


-1.1 


8.34±0.15 


0.5 


18±6 


11.02±0.08 


-1.1 


8.32±0.11 


0.7 


16±5 


11.09±0.15 


-1.1 


8.33±0.10 


0.9 


12±4 


ll.08i0.ll 


-1.1 


8.17±0.18 







Red sequence 






0.0 


37±4 


10.81±0.10 


-0.7 ±0.07 


8.33±0.10 


0.3 


17±7 


10.97±0.09 




-0.7 


8.17±0.18 


0.5 


15±5 


10.95±0.10 




-0.7 


8.08±0.11 


0.7 


11±4 


11.06±0.18 




-0.7 


8.05±0.10 


0.9 


9±3 


11.01±0.08 




-0.7 


7.89±0.18 



Blue cloud 



0.0 


9.4±1.0 


10.80±0.10 


-1.45 ±0.03 


7.98±0.10 


0.3 


7.5±3.6 


10.86±0.19 


-1.45 


7.94±0.15 


0.5 


7.2±2.2 


10.93±0.12 


-1.45 


8.00±0.11 


0.7 


6.4±2.2 


11.07±0.07 


-1.45 


8.08±0.10 


0.9 


5.3±1.7 


11.00±0.25 


-1.45 


7.93±0.10 



Focusing instead on the red sequence (diamonds and 
solid lines), there is rapid evolution in their stellar mass 
function. Assuming a constant faint-end slope (which we 
cannot meaningfully constrain in this work^), the formal 
fits prefer evolution in </>* (such that (j)* increases towards 
the present day) without significant change in M* . It is 
worth noting that this result is insensitive to the detailed 
choice of color-mass cut used to separate the red and blue 
populations — as long as the cut is somewhere in the 'gap' 
between the red and blue populations. 

4.5. The evolution of stellar mass densities 

The average stellar mass density p^ = J ^(m)dm is cal- 
culated by integrating the Schechter functions of Table 
n between m = and oo. The mass density p* for all 
galaxies is calculated for all three COMBO-17 fields; the 
larger of i) the field-to-field scatter, divided by -\/2, or ii) 
the predicted variance, following Somerville et al. (2004) 
given a number density and redshift, is adopted as the 
cosmic variance error bar. The results in Fig. ^| show, 
that the total stellar mass has increased substantially be- 
tween 2; ~ 1 and the present day. Interestingly, much of 
the evidence for this growth comes from the comparison 
to the local stellar mass density^. Importantly, this stellar 



* Kodama et al. (2004) and de Lucia et al. (2004) present 
deeper data which show some evidence for a deficit of faint 
red-sequence galaxies at 2; ~ 1 compared to the present-day 
red-sequence galaxy stellar mass function. 

^ In this context, it is worth noting that the local total stel- 
lar mass densities, and indeed those split by color, calculated 
by a variety of different groups are very similar with < 0.1 dex 
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Fig. 10. The integrated stellar mass density as a function 
of redshift. In the upper two panels the total mass den- 
sity for all galaxies (filled circles) is compared with those 
for red-sequence galaxies (diamonds), and for blue cloud 
galaxies (asterisks) seperately. The lower panel shows the 
fraction of mass in red sequence galaxies as a function of 
redshift. In all cases, mass functions are integrated down 
to zero mass and error bars come from field-to-field vari- 
ation divided by ^/2. The z = datapoint is taken from 
Bell et al. (2003). 

mass growth is dominated by the growth of the total stel- 
lar mass in red-sequence galaxies, while the total stellar 
mass in blue cloud galaxies is roughly constant over the 
last 8 Gyr. 

5. Discussion 

The COMBO-17 data and analysis constitute one of the 
first opportunities to combine large samples of galaxies to 
z^ 1, rather precise redshifts (average cr^/(l ± z) ~ 0.02), 
and accurate rest-frame colors (U-V)rest luminosities and 
stellar mass estimates. In particular, its rest-frame col- 
ors are accurate enough to allow for the first time the 
exploration of the mass growth of red and blue galaxies 
separately. 

5.1. SED-based Stellar Mass Estimates 

Through extensive SED matching, based on the PEGASE 
models (Fioc and Rocca-Volmerange 1997), we have illus- 

scatter if IMFs and cosmologies are accounted for, giving con- 
fidence in the robustness of the locally-measured stellar mass 
density estimates (e.g.. Cole et al. 2001; Kochanek et al. 2001; 
Bell et al. 2003; Baldry et al. 2004; Panter et al. 2004). 
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Fig. 11. A comparison of stellar and dynamical masses. 
Accurate stellar and dynamical masses (asterisks) are 
taken from Cappellari et al. (2005) for a sample of lo- 
cal early-type galaxies. Dynamical masses (van der Wei et 
al. 2005) for 6 COMBO-17 galaxies at 0.6 < z < 1.05 with 
robust stellar masses are shown as diamonds. 



trated the strength and limitations of stellar mass esti- 
mates based on this approach. We have shown that for 
a given IMF (e.g. Kroupa et al. 1993), the present mod- 
elling indicates that SED-based M / L estimates are better 
than a factor of two for most individual galaxies, with sys- 
tematic sample uncertainties probably < 30% (c./. Fig.^] 
and Section 3.5, see also Drory et al. 2004b). We stress, 
however, that the masses derived here simply scale for dif- 
ferent IMFs (below 2Mq); e.g. for a Salpeter (1955) IMF 
all masses would be higher by a factor of 1.8. 

Note that the analysis of the redshift dependence of 
stellar masses does not require an accurate absolute cal- 
ibration of the stellar mass scale; rather, relative consis- 
tency between the stellar mass scale at 2 ~ 1 and 2 ~ 
is sufficient. Figure 11 shows the comparison between a 
sample of local early- type galaxies with accurately-derived 
stellar and dynamical masses (Cappellari et al. 2005; a 
Kroupa 2001 IMF is adopted in that work, which has a 
very similar mass scale to the Kroupa et al. 1993 IMF) and 
a sample of 6 distant early-type galaxies in the Chandra 
Deep Field South with robust dynamical mass estimates 
from deep spectroscopy (van der Wei et al. 2005) and stel- 
lar mass estimates from this paper^°. There is a clear 



A substantial fraction of the van der Wei et al. sample were 
at or below the COMBO-17 tur < 24 limit and had poorly- 
estimated or undefined redshift. 



relationship between stellar and dynamical mass^^, and 
this relationship seems independent of redshift at least 
out to 2; '--^ 1. This gives us some confidence that stellar 
population masses are being estimated consistently for lo- 
cal and distant samples, thus evolutionary trends between 
COMBO-17 and local samples should be relatively robust. 
Our study illustrates quite dramatically how galaxy 
samples, even when selected at red observed bands, are 
very far from a stellar mass limited sample. The difference 
in stellar mass between the bluest and reddest galaxies at 
a given flux limit is a factor 8 at 2 ~ 0.3 and more than 30 
at 2 ~ 1 (where the sample-selecting R-band corresponds 
to Arost ~ 360 nm). 

5.2. The relationship between color and stellar mass 

Our most conspicuous astrophysical result is that at z < 1, 
the majority of the most massive galaxies (M* > IO^^Mq) 
are on the red sequence (see Figs. 7 and 9). This has 
qualitatively been clear from color-luminosity diagrams 
(e.g., Bell et al. 2004), but Figs. 7 and 9 demonstrate 
that it holds true for the top decade in stellar mass. It is 
not clear, at this stage, what this result implies in terms 
of galaxy evolution. The stellar mass functions (Fig. 9) 
show a very modest increase in the number of red massive 
(M, > lO^^Mo) galaxies and a modest decrease in the 
number of blue massive ones from 2 = 1 to the present 
day. Yet, it is important to bear in mind that uncertainties 
in stellar masses and photometric redshifts have a partic- 
ularly strong effect on the high-mass end of stellar mass 
functions; therefore, we argue that it is not clear whether 
the massive galaxy population was largely in place and 
non-star- forming at 2 ~ 1 (e.g. Saracco et al. 2005), or 
whether there has been some growth of massive non-star- 
forming galaxies through processes such as dry mergers 
(e.g. BeU et al. 2005b). 

Over the last decade, there has been extensive evidence 
for what has become known as 'downsizing': while mas- 
sive galaxies form stars at high redshift (e.g., Chapman et 
al. 2003; Daddi et al. 2005), only low-mass galaxies form 
significant numbers of stars at the present (Cowie et al. 
1996; Heavens et al. 2004; Juneau et al. 2004; Bauer et 
al. 2005; BeU et al. 2005a; Tanaka et al. 2005; Kelm et al. 
2005). The color-split mass functions presented in this pa- 
per allow us to explore some aspects of this effect. Figure 
9, for instance, shows the evolution of the mass at which 
the red sequence {— non-star-forming) and blue cloud (= 
star forming) stellar mass functions cross^^ very similar to 
that adopted by Kauffmann et al. (2003), who explore the 

^^ Note that dynamical mass should always be in excess of 
stellar mass (allowing the existence of dark matter in the inner 
parts of galaxies); clearly a Kroupa et al. (1993), Kroupa (2001) 
or Chabrier (2003) IMF satisfies these constraints whereas the 
Salpeter (1955) IMF (which produces stellar masses ~ 0.25 dex 
larger) would violate this constraint (see Cappellari et al. 2005; 
BeU & de Jong 2001). 

^^ This definition is identical to Bundy et al.'s (2006) transi- 
tion mass Mtr, and is 
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distribution of 4000 A break strengths as a function of stel- 
lar mass.: above this mass the majority of galaxies is red 
and below it the majority of galaxies is blue. We find that 
this characteristic mass is ^ 10^"'*Mq at z = 1 (in agree- 
ment with estimates presented in Bundy et al. 2006) and 
decreases by a factor of 5 or more to ~ IO^^Mq at z ~ 0. 
Yet, the physical significance of this 'transition mass' is 
not obvious as the stellar mass function of blue galaxies 
does not change significantly in the epoch < z < 1, 
suggesting that the maximum mass that galaxies reach 
while forming stars remains roughly the same. The appar- 
ent evolution of the 'transition mass' is driven entirely by 
the increasing prominence of the non-star-forming red se- 
quence galaxies. This emphasizes the need for care when 
discussing 'downsizing', and the utility of large samples of 
galaxies with robust estimates of stellar masses for study- 
ing the physics driving the order-of-magnitude decrease in 
cosmic star formation rate since z ~ 1 (e.g., Le Floc'h et 
al. 2005). 

5.3. The build-up of stellar mass 

It is interesting to compare the COMBO- 17-derived stellar 
mass density evolution since z ~ 1 with other work that 
has studied the build-up of the galaxy population since 
z ~ 4. The results are shown in Fig. 12. The top panel 
shows the total stellar mass in galaxies as a function of 
redshift (integrated to zero mass) for a variety of studies. 
The stellar masses have been adjusted to be consistent 
with a Kroupa et al. (1993), Kroupa (2001) or Chabrier 
(2003) IMF; all three IMFs have a stellar M/L at a given 
optical color that are 0.25 dex lower than for a Salpeter 
(1955) IMF, and all three IMFs are consistent with dy- 
namical constraints (Bell & de Jong 2001; Cappellari et al. 
2005), while a Salpeter IMF over-predicts the dynamical 
masses of many galaxies. Despite the inevitable uncertain- 
ties associated with estimating stellar masses, and cosmic 
variance. Fig. 12 illustrates that the different studies have 
come to a reasonably consistent picture of the growth of 
stellar mass with cosmic time. 

This growth of stellar mass should be consistent with 
the integral of current measurements of the cosmic SFR. 
The lower panel of Fig. 12 shows the evolution of the 
cosmic SFR (adapted from Hopkins 2004). We have cho- 
sen not to convert the cosmic SFR measurements from a 
Salpeter (1955) IMF to a Kroupa et al. (1993) IMF in this 
work: as the power-law index of a Kroupa et al. (1993) 
IMF is substantially steeper than those of a Salpeter, 
Chabrier or Kroupa (2001) IMF for high-mass stars (which 
dominate the SFR measurements, but do not significantly 
affect broad-band SED-derived stellar masses), the use of 
a Kroupa et al. (1993) stellar IMF in what follows would 
lead to a factor of three over-production of stars. Instead, 
we adopt a Chabrier (2003) IMF (which has the same 
high-mass IMF slope as a Salpeter IMF and very similar 
stellar M/Ls to a Kroupa et al. (1993) or Kroupa (2001) 
IMF), and ask the question whether the integral of the 
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Fig. 12. Top: the growth of stellar mass density with cos- 
mic time, as estimated from a number of surveys. Bottom: 
the evolution of the cosmic star formation rate, adapted 
from Hopkins (2004). In both panels, the stellar IMF has 
been adjusted to be consistent with Chabrier (2003); a 
Kroupa et al. (1993) IMF has a stellar M/L consistent 
with such an IMF but would produce star formation rates 
a factor of three higher than the Chabrier (2003) IMF. The 
lines show different assumptions for the evolution of the 
cosmic star formation history and the associated growth 
of cosmic stellar mass: the dotted line shows the fit opti- 
mized to reproduce the growth of stellar mass, the dashed 
line fits the star formation rate better, and the solid line 
is a reasonable compromise between the two. 



cosmic star formation rate is consistent with the observed 
growth of stellar mass. We show three cosmic star for- 
mation histories and their integrals (accounting fully for 
gas recycling). The dashed line is a fit optimized to repro- 
duce the evolution of cosmic star formation rate: this tends 
to produce a little too much stellar mass at all redshifts. 
The dotted line provides a much better fit to the build-up 
of stellar mass, but appears to undershoot somewhat the 
star formation rate estimates at z > 1. The solid line is a 
reasonable compromise between matching the cosmic star 
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formation rate and build-up of stellar mass. Despite these 
interesting tensions at the 30-50 % level — which may 
indicate deficiencies in stellar mass or SFR estimates {e.g. 
due to a globally incorrect IMF and/or one which varies 
from galaxy to galaxy), or may betray inappropriate ex- 
trapolation to total densities from the observed galaxies — 
it is nonetheless fair to say that the match between cosmic 
SFR and the build-up of stellar mass is rather good. 

The relatively precise rest- frame colors of COMBO- 17 
allowed us to separate out the evolution of (yO*)(z) in red 
and blue galaxies. While a similar split has not yet been 
carried out for estimates of cosmic star formation rate, it is 
interesting to note that study of the UV (Wolf et al. 2005) 
and UV-I-IR (Bell ct al. 2005a) luminosity of galaxies has 
demonstrated that the bulk of star formation takes place 
in spiral galaxies residing in the blue cloud. Thus, one 
would naively expect that the growth of stellar mass at z < 
1 would be dominated by a growth in the total stellar mass 
in blue galaxies. Yet, Fig. 10 shows that the total stellar 
mass growth takes place mainly in the red sequence. Thus, 
stars form in blue galaxies, but a substantial fraction of 
the star-forming galaxy population must globally truncate 
their star formation and fade and redden onto the red 
sequence at z < 1. A discussion of this point is deferred 
to a future paper. 



6. Summary and Conclusions 

We have presented estimates of the stellar masses, derived 
from SED-fitting of the COMBO-17 flux-points using 
PEGASE population synthesis models, for 25 000 galax- 
ies to z = 1, with known redshifts, rest-frame colors and 
luminosities. From this data set we were able to show the 
following. 

— Since z ~ 1 the majority of galaxies with M > 10^-^ Mq 
have been on the red-sequence, i.e., have essentially 
completed the build-up of their stellar mass. 

— As a consequence, the stellar mass function of red and 
blue galaxies differ strongly from each other at all red- 
shifts from to 1. Red galaxies always dominate the 
massive end, while blue galaxies dominate the faint 
end. The mass at which the mass functions cross de- 
creases from M* ~ lO^^'-^AfQ at z = 1 to ~ lO^^Mg at 
the present day. This evolution in the 'transition mass', 
above which most galaxies are not forming stars, is 
driven primarily by the increasing prominence of the 
non-star-forming galaxy population at later times. 

— The total stellar mass in galaxies has roughly doubled 
in the 8 Gyr since z = 1. Assuming a Chabrier (2003) 
IMF, the stellar mass evolution and star formation rate 
evolution are reasonably consistent: a Kroupa et al. 
(1993) IMF (which has a stellar M/L for a given SED 
which is very similar to those of a Chabrier IMF) has 
a relatively low number of high-mass stars and no con- 
sistent interpretation of the cosmic star formation rate 
and growth of stellar mass can be found. 



— Somewhat counter-intuitively, the stellar mass func- 
tion of blue, star-forming galaxies has stayed more or 
less constant since z = 1; the mass in "dead" red- 
sequence galaxies has almost tripled over that period. 
While at z = 1 half the stars lived in red sequence 
galaxies and half in star-forming galaxies, 2/3 of all 
stars at the present epoch live in red-sequence galax- 
ies. 

— Even with 25 000 galaxies in three disjoint fields of 
30' X 30' each, field-to-field variations limit the pre- 
cision of (p*)(z) estimates to ~ 30%. 

There are a number of obvious steps which need to be 
taken to improve the accuracy and extend the scope of 
the work presented in this paper. Larger areas, more pre- 
cise redshifts (particularly spectroscopic), and wider wave- 
length coverage will give increased precision on the form 
of the mass function (especially at the high-mass end; crit- 
ical for exploring the infiuence of galaxy mergers on the 
stellar mass function) and increased redshift range. Larger 
samples of dynamical masses and improved templates will 
increase the accuracy (i.e., reduce systematic error) in the 
mass functions. Finally, accurate determinations of star 
formation rates from Spitzer/Herschel/ALMA will allow 
star formation and the build-up of stellar mass to be inter- 
compared at a much more detailed and informative level 
than is currently possible, further elucidating the role of 
in situ star formation vs. galaxy mergers in driving the 
mass evolution of the galaxy population. 
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